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A free Fermi gas has, famously, a superconducting susceptibility that diverges logarithmically at 
zero temperature. In this paper we ask whether this is still true for a Fermi liquid and find that 
the answer is that it does not. From the perspective of the renormalization group for interacting 
fermions, the question arises because a repulsive interaction in the Cooper channel is a marginally 
irrelevant operator at the Fermi liquid fixed point and thus is also expected to infect various physical 
quantities with logarithms. Somewhat surprisingly, at least from the renormalization group view- 
point, the result for the superconducting susceptibility is that two logarithms are not better than 
one. In the course of this investigation we derive a Callan-Symanzik equation for the repulsive Fermi 
liquid using the momentum-shell renormalization group, and use it to compute the long-wavelength 
behavior of the superconducting correlation function in the emergent low-energy theory. We expect 
this technique to be of broader interest. 



I. INTRODUCTION 

Thinking within the framework of the renormaliza- 
tion group (RG) is always an insightful way to approach 
the low energy, long wavelength behavior of physical 
systems, but it docs not always provide a significant 
computational advantage over more straightforward, per- 
turbative, methods. The exceptions are cases that in- 
volve marginal and near marginal couplings where mak- 
ing sense of pcrturbativc divergences is far easier within 
the framework of the RG. Indeed, the historical develop- 
ment of RG methods had mostly to do with exactly these 
cases leading up to the discovery of the epsilon expansion 
by Wilson and Fisher— 

A subclass of these canonical RG applications con- 
sists of problems with marginally relevant or irrelevant 
couplings. Celebrated examples of the former are the 
Kondo problem^ and QCD^ while the textbook exam- 
ple of the latter is the critical behavior of ferromagnetic 
models (equivalently, of vector models in the language of 
field theory) in four dimensions^. In all of these cases, 
non-trivial logarithms appear in physical quantities ei- 
ther in the ultraviolet (marginally relevant) or the in- 
frared (marginally irrelevant), and reflect the extremely 
slow variation of the coupling at issue. An instructive ex- 
ample, with close parallels to our concerns in this paper, 
is the critical behavior of the specific heat in four di- 
mensional O(N) ferromagnets/vector models. Here the 
non-interacting gaussian theory exhibits a logarithmic di- 
vergence 

T — T 

C log \t\ with t = . 

T c 

As the quartic coupling is marginally irrelevant, i.e., it 
dies logarithmically in the infrared and the theory flows 
to the non-interacting limit, one may then suppose that 
the divergent susceptibility of free-field theory will re- 
main intact even upon inclusion of the interactions. How- 
ever, as is long established, interactions lead to the more 



complex logarithmic dependence 

C~[-log|t|]*3 

which diverges as (— log l^l) 1 ^ 3 in the Ising case (N = 1) 
but vanishes as (— log when N — > oo. Observe that 
the logarithms appear generically with fractional powers 
and a variable sign. This complexity is not easily unrav- 
eled pcrturbativcly and instead it is much better to resort 
to the machinery of the RG, specifically to the derivation 
and solution of the RG (differential) equations obeyed 
by the correlation functions, which we shall simply call 
Callan-Symanzik equations in this paper. 

Let us turn now to the problem considered in this pa- 
per. The RG treatment of interacting fermions identifies 
the Landau Fermi liquid as a fixed point of a momentum 
shell RG characterized by set of exactly marginal cou- 
plings consisting of the Landau F function. In addition 
however, there is a marginally irrelevant coupling which 
is the repulsive BCS (Cooper channel) coupling^ Absent 
this marginal flow the fermions have a logarithmically di- 
vergent susceptibility to superconductivity. We wish to 
ask what happens to this divergence when the marginally 
irrelevant flow is taken into account. We show here us- 
ing the RG that this divergence goes away and the uni- 
form, zero frequency superconducting susccptibilty of the 
Fermi liquid is therefore finite. 

It is useful at this point to clarify what we mean by 
the Fermi liquid. As the reader is no doubt aware, the 
actual RG flow for an interacting fermion system with 
typical repulsive interactions leads inevitably to a su- 
perconducting instability via the Kohn-Luttingcr effect 
wherein screening produces an effective interaction that 
is attractive in a higher angular momentum channel, as 
discussed for example in Ref 0- What we have in mind 
therefore is the RG flow for a system where the dom- 
inant bare couplings in the BCS channel are repulsive 
and renormalize to smaller values for a large range of 
scales (temperature or energy) while the growing attrac- 
tive couplings arc still small — operationally this is what 
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one means by a Fermi liquid. We can, however, formalize 
this understanding by working with Hamiltonians which 
contain only the Landau couplings and the reduced BCS 
couplings that are left in a naive application of the RG to 
the interacting fermion problem. In this approximation, 
the Kohn-Luttinger scale has truly vanished but we are 
still left with a marginally irrelevant operator about the 
Fermi liquid Hamiltonian whose flow can be studied. 

Our central result is possibly not new — certainly it is 
almost present in the large lore on superconductivity and 
superconducting fluctuations^ and we would be delighted 
to hear from readers who can point us to an explicit, rel- 
evant citation. It does not, however, appear to be widely 
known and at first blush comes as a surprise, as con- 
densed matter theorists are conditioned to think of the 
Fermi liquid as exhibiting correlations morally identical 
to that of the Fermi gas up to the effects of the Landau 
parameters. Ex post facto the intuition we would offer 
is that the repulsive interaction goes away just slowly 
enough at long distances that the usual buildup of the di- 
vergence of the superconducting susceptibility is undone. 
Possibly there is a deformation of our present problem 
where this exact cancellation can be modified to yield a 
more complex residual of the kind cited for four dimen- 
sional ferromagnetism above. 

In any event, we view our work as a contribution to 
the RG analysis of interacting fcrmions and as such we 
trust readers will find it interesting as well. In the pro- 
cess, we show how to derive Callan-Symanzik equations 
to one loop for interacting fcrmions for composite opera- 



tors made from the fundamental Fermi fields — which may 
be of interest to readers even beyond the specifics of our 
application to the computation of the superconducting 
susceptibility. 

In the following we begin with a quick summary, in Sec- 
tion II, of the RG for interacting fermions that leads to 
the Fermi liquid fixed point (s) and the flow in its vicinity. 
Next we derive the relevant Callan-Symanzik equations 
(Section III), solve them (Section IV) and end with a 
brief discussion (Sections V and VI) and an Appendix 
that contains some technical details. Note that while 
our results certainly apply in both two and three dimen- 
sions, to avoid unnecessarily complications we work with 
the pedagogically simpler case of d = 2 throughout this 
paper. 



II. REVIEW OF RESULTS FROM THE RG 

We begin the technical part of our discussion with a 
summary of results from the renormalization group ap- 
proach as applied to interacting fermion systems. We 
shall provide a telegraphic review, referring the reader 
interested in further details to more pedagogical discus- 
sions, such as Ref. 0- Readers familiar with the technol- 
ogy and results can skip ahead to Section III. 

We focus on the following following action, written for 
a system of electrons with a circular Fermi surface in 
d = 2 and spin directions a =1, ! : 



= So + Si 



where ^> M (i) = if> '^(uidiki) is the Fermion/ Grassman field, 
k = |K| — Kp is the radial component of momentum 
measured relative to the Fermi momentum Kp , vp is the 
Fermi velocity, A << Kp is the ultraviolet cutoff. 
In the second term Si, the measure is 

1 *' K 4 =Ki+K 2 -K 3 

(2) 

The renormalization group transformation involves 
three steps: (i) integrating out all momenta between A/s 
and A, and correcting terms in the action as needed in 
the process; (ii) rescaling frequencies and the momenta as 
per [ui, k) — ► s (w, k ) so that the cutoff in k is once again 
at ±A; and finally (iii) rescaling fields ip — > Z^ip to keep 
the free- field action So invariant. Thus So is a fixed point 



t (u0k) + ^2 )V ^(l, 2,3,4)^(1)^(2)^(3)^(4) 

(1) 

I 

of this RG and the possible interactions can be classified 
as relevant, irrelevant, or marginal with respect to this 
transformation and fixed point. 

At tree level, it is easily shown that interactions with 
six or more fields are irrelevant, as is the a; or A: depen- 
dence of quartic couplings. While this is very much like 
<\>\, the 4 theory in four dimensions, there the coupling is 
•just a number, u(0, 0, 0, 0), describing collision of particles 
at zero external momentum, here the quartic couplings 
can depend on the angles on the Fermi circle. Given that 
momentum is conserved, we need to pick only three of 
the angles independently, say 8i,9 2 and 6 3 . However, 
the fact that the momenta come not from the plane but 
a very thin annulus leads to additional constraints 7 . 

Consider the left half of Figure Q] where all momenta 
lie on the Fermi circle. Given 1 and 2 , it is clear that 
3 and 4 have to equal them pairwise. That is 0\ = 63 
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and 62 = 84 or the exchanged version. Now consider 
the generic case of momenta that lie in a very thin shell 
|fc| < A, rather than right on the Fermi circle. It is not 
surprising that now 8\ ~ 83 and 82 — &4 with deviations 
of order AjKp. We ignore the dependence of u on such 
tiny angular differences when A/Kf — > 0, and we ignore 
any k -dependence (since it is irrelevant). Thus we define 
u with all k — and 8\ = 83 and 82 = 84. 

u(8 1 , 82,83,84) = u(8i,8 2 ,-8 1 ,-9 2 )=u(8 1 ,8 2 ) 

= u(8 1 -92) = F(8). (3) 

It is important to bear in mind that even though we 
evaluate u at k — for angles {8\ = 83) and (82 =84), 
we do not imply that only forward scattering is allowed 
or that all momenta lie on the Fermi surface: rather we 
allow all momenta in the measure defined by Eqn. [2] 
but ignore the dependence of u on k and the tiny differ- 
ences (81 — 83) and (6*2 — 84). (It is like saying that in 
</>|, u = u(0, 0,0,0) does not mean the external legs are 
limited to zero, only that u is the same for all values of ex- 
ternal momenta.) Indeed a small amount of non- forward 
scattering is not only allowed at any nonzero cut-off, it is 
essential to produce a nonzero compressibility in the "g" 
limit. 

In summary, in studying the flow of the four point cou- 
pling, we will choose the legs to be on the Fermi surface 
(since k dependence is irrelevant), i.e, we study the flow 
off. 

Kinematics allows one more coupling function besides 
F. Consider the right half of Figure [T] with all momenta 
on the Fermi circle, but with 9± = —9 2 so that the incom- 
ing momenta add up to zero. Now the outgoing pair of 
momenta can point in any direction, as long as they add 
up to zero, which lets them point in any pair of mutually 
opposite directions. This leads to a function V in the 
BCS channel: 





FIG. 1: Kinematics of the couplings F and V. 



in V can be accounted for by expanding in different an- 
gular momentum channels, and deriving flow equations 
for each channel. Each channel will have its own Callan- 
Symanzik equation, and the remainder of the arguments 
follow from this. On the other hand, it should be clear 
from the preceding discussion that restricting V to be 
repulsive is necessary, since attraction leads to the BCS 
instability, and of course the Fermi liquid fixed point can- 
not capture the physics of superconductivity 
Thus. 
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^ M (1)^(2)^(3)^(4)6(A-|P|) 



(5) 

where the additional condition 0(A — |P|) ensures that 
the total momentum of the Cooper pair is less than A 
which itself is assumed to obey A << A. This will not be 
a limitation in our study of superconductivity since only 
the interaction at P ~ will come into play. 
We will use the standard parametrization of s: 



u(6 u 82,63,84) = u(0 1 ,-6 1 ,6a,-0 a )=u(0i,6 a ) 

= u(9 1 -6 3 ) = V(8) (4) 

Once again, even though u will be assumed to depend on 
just the directions 8\ and 83 (via their difference), for any 
A > 0, scattering between states of the total momentum 
P ~ A will be kinematically permitted. 

The fate of the couplings F and V, marginal at tree 
level, is determined by one loop diagrams. The one-loop 
analysis of Ref shows that F is strictly marginal while 
V flows. It is also shown there that the one-loop flow is 
exact in the limit A/Kp 0. 

In the remainder of this paper, we will assume F = 
and that V is repulsive and angle independent. While 
this simplifies the discussion, the first and last restric- 
tions are inessential. Inclusion of F will not affect the 
physics of the Cooper channel for kinematic reasons, and 
while the angular dependence of V can affect details of 
various calculations, it cannot alter the fact that singu- 
lar logarithms arise in the same. As long as the system 
remains rotationally invariant, any angular dependence 



_ . , ds dA 
s = e so that da = — = — — . (6) 
s A 

We will typically choose to rcscale by an amount close to 
unity, in which case 



da = — = ds since s — > 1. 
s 



(7) 



Note that the cutoff is being changed infinitcsimally at 
each RG step; as the RG transformations arc infinitesi- 
mal, this permits us to derive differential equations that 
capture the content of the flows. 
It is shown in Ref. that 



0, \ - d 9 2 
P{9) = t = ~ag 
da 



(8) 



where a is a positive constant. The solution to this equa- 
tion is 



9o 



1 + aag 



(9) 
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For repulsive go we see a logarithmically vanishing g(a) ~ 
1/ct. One may expect that such a marginally irrelevant 
coupling cannot affect the superconducting susceptibil- 
ity which diverges logarithmically in the noninteracting 
limit. The rest of this paper aims to show why this is not 
the case and what finally happens. 



III. DERIVATION OF THE 
CALLAN-SYMANZIK EQUATIONS 

So far, we have reviewed how different coupling func- 
tions in the Hamiltonian flow under the RG. This allowed 
us to classify relevant and irrelevant perturbations, and 
determine that the repulsive Cooper channel interaction, 
being marginally irrelevant, does not lead to an instabil- 
ity: the system flows. In this section, we discuss how to 
derive a set of differential equations that describe how the 
correlation functions evolve under the RG flow. First, we 
discuss how to derive the equations for correlation func- 
tions without composite operators, and then discuss the 
procedure for the Cooper pair operator. 



Callan-Symanzik Equations for fundamental 
Fermi field correlators 



We begin with the correlation function-^, defined with 
an action S and cutoff A: 



[D-ijjDip] j 



-Sa[4>M 



^aSk^ek 



[D-ipDip] A e 



-Sa[1>,4>] 



(10) 



where [DvfjDip] A indicates that the functional integral is 
over modes of all momenta upto A. Here it is understood 
that the momenta k are well within the cutoff. Note also 
that the momenta and frequencies corresponding to tp 
and ip are different; in fact, the result will have a delta 
function that forces these to be the same, but since the 
delta function is dimensionful, it affects the scaling, so we 
maintain different momenta explicitly through the calcu- 
lation. Integrating out a shell dA of momenta between A 
and A + dA = A - |gJA| = A/s = A/(l + da), we haveZ. 



(ipsjSk^e^sA 



Qek^Sk 



-s="[v- < ^ < ] 



(11) 

where S^/f is an action for a theory with cutoff A/s, 
with the appropriate corrections to parameters, and the 
< superscript denotes 'slow' modes that lie within the 
reduced cutoff. In order to obtain a theory with the 
same cutoff, we rescale momenta in the second equation. 
We will write ip'(uj'6'k') = Z^ 1/2 ^<{u' / s,6,k' / s), where 
< |fc'| < A. Further, we assume that action may be 



written as S e A f / s [^<, $<] = S A [tp', f) + SS A [tp', tp'}, where 
SS is O(da). We obtain 



ui6k/S,A 



[D^D^} A e- s ^ ss ^ Ss - k ^ su]9sk 



[D^£>V'] A e- 



-Sa+<5Sa 



Z^{lp sG j t g^ki>suj,6,sk)s+SS,A 



(12) 



where in the second step we dropped the primes as the 
fields are integrated over. 

Since SS depends on the parameter s of the flow, we 
can write S + SS = S s ; then, we have 

{i>QGk^ek)s,A = Z 4 ,(ij sQgsk -ifi SUJ ,0, sk ) SsA (13) 

The 1 left side- of this equation is manifestly independent 
of s. Therefore, differentiating both sides with respect to 
s, 

and dividing through by Z^, 



J— -(^sa,6,ski J suj.d.sk)s s ,A + -r(' l PsLD,0,sk' l Psuj,9,sk)s 3 ,A - 

ds ds 
The second term can be rewritten aa^ 



ds 



(i> S Q,e,ski'su 1 ,e,sk)s s ,A 
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duo dk 



dg_d_ 

ds dg 



su,e,sk/S s ,A 



(14) 



where we note that the only change in the action to one 
loop order is in the coupling constants. Collecting terms, 
and taking the s — > 1 limit, we arrive at the Callan- 
Symanzik equation for the two-point correlation function, 



d_ 

duj 



8 



d 



k dk + ^dg- +3 - 2 ^ 



where f3{g) = % 



{^QSk^eklsA = 
(15) 

4 s is the beta function, and we 

da 1 



define the anomalous dimension 3 — = — 4— 

Next, we observe that the Green's function is re- 
lated to the correlation function as (ipQgkiJuek) S,A = 
G(u), 9, k)S u ^6 e §S k f, . Since the 5 functions each have 
dimension —1, it is easily verified that the Green's func- 
tion satisfies the relation 



G{uj,9,k) = (16) 



Finally, we rewrite this as an equation for the ampu- 

G-h 



tated Green's function, T^ 2 ' ' 1 



d , d n . . 8 



dg 



r^ o) ( w ,0,fc) = o 



(17) 
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We see that for the free theory, where 



.(2,0) 
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0, 



iui — vpk satisfies the above relation, as expected. 



B. The Callan-Symanzik Equation for the Cooper 
Pair Operator 

The central objects of this paper are the Cooper pair 
operator and its two-point correlation function. Wc 
define the composite operator that creates an s-wave 
Cooper pair at frequency Q, and momentum P << A 
by 



On: 



J -co 2^. 



r A dk 

—A 2tt 



r 2ir 



d9 - 

2^< 



(18) 



where (k, 9) refer to K. (Since P « A, if K lies in the 
shell so will K + ^.) The operator O is likewise made of 
two ip's. 

Consider the expectation value 



(19) 

In the first stage of the RG transformation, we must 
integrate out all fields at momenta between A/s and A, 
so that the remaining action only has terms that are at 
momenta < A/s. This is a complicated process for the 
composite operator since it contains terms that must be 
integrated out. We can write the composite operator as 



/ 


~[ DtpDip 

. A 




/ 


_ A 


e -s A li>,i>] 



o 



= o, 



n.p 



a <,> 



>,< 

W Q,P 



>,> 



(20) 



where the < , > denote whether the two fields entering 
the operator are below or above the cutoff respectively. 
The operator O^'p is the descendant of the composite 

operator in the theory with cutoff A/s; Oqp and p 
are mixed terms, whose > field must be integrated over; 
and Oqp is composite operator made up entirely of fast 
modes that will be integrated out. In integrating over 
the fast modes we may functionally average over the free 
field action, since the deviations from the free action will 
only produces terms of higher order in g. 

When we take the product of composite operators, 
we note that only terms with even numbers of 'fast' 
(>) modes will survive the functional average over fast 
modes. Of these, one is the term O^'pO^p , which 
is the descendant of the two-point function in the 
lowered-cutoff theory; then there are the mixed termsi^, 
O^'pOn'p and O^'pO^'p*; and finally there is the piece 

made up entirely of fast modes, OnaOpp, The first 



term is analogous to i/> < i/} < in our earlier discussion of 
the fermion Green's function. The remaining three pieces 
have no analog for non-composite operators, and we must 
determine how they alter the RG flow. 

They clearly form additive corrections to what comes 
from the descendant term. Just for computing this ad- 
ditive term we simplify things by invoking a familiar re- 
sult from computing the T = Cooper bubble explicitly: 
namely, that the correlations depend on (P, 0) only via 
the combination fl = \/f2 2 + P 2 . So we will simplify the 
following discussion by choosing P = 0. 

Now there are no mixed terms since | K + y | = | — K + 
1^1 (when P = 0) and both fields in Oq and O n are 
cither above the cut-off or below it. 

This leaves only the term consisting entirely of fast 
modes. It remains to integrate over K, uj, i.e. we have (on 
converting to the Fermi-surface coordinates, performing 
the trivial angular integral, and recalling that all stray 
factors of Kp from the measure are absorbed into the 
definitions of the fields) 



(O 



>,> ro >,> N 



dk 

-> 2 x / — 



A/s 



2- 



_^ 2?r uj 2 + vlk 2 



1 



. 1 - - 

2nvp \ s 



(21) 



As was mentioned earlier, this is a term of order 
— dA/A=l — l/s = dcr;we show below that this leads to 
an inhomogeneous term in the C-S equation for the two- 
point function. As an aside, we point out (but do not 
prove) that for any higher correlation function of com- 
posite operators beyond the two-point function, and for 
any other correlation functions with external legs and 
composite operator insertions, this contribution is of a 
higher order in da, and will not lead to inhomogeneous 
terms. For aficionados of the RG for the </f> 4 theory, this 
is the fermionic analog of the statement that the spe- 
cific heat is the only operator that is not multiplicatively 
renormalizable, the signature of which is the appearance 
of inhomogeneous terms in the Callan-Symanzik equation 
for the specific heat£. 

Finally, we turn to the more familiar term, which is 
the 'descendant' operator in the cutoff theory. We still 
need to determine how this behaves under the next two 
RG steps. It is sufficient to consider how the compos- 
ite operator gets rescaled; the two-point function scaling 
follows immediately. Before implementing the RG trans- 
formation, we should rewrite the composite operator in 
terms of the momentum relative to the Fermi momentum. 
Using the fact that K = (K F + fc)K, where k e [—A, A] 
at any stage of the RG, on expanding to linear order 
in the momentum P, we find the displacement from 
the Fermi surface is given by k± = |±K + ^ — Kp w 
k ± PcosOkp, while the angular coordinates are given 
by 9 + = 9 k - g|U sin 6 kP ,6- = ir + 9 k + ^-sm9 kP . 
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Therefore, 



It follows that the descendant operator is 



O 



fe,OJ 



ifo w + Q, k + Pcos 6» fcg , 0fc 



P 
2iO 



■ sin 6» fc p 



P 



x ^4. 1 -w,k- Pcos6 k p,w + 6 k + 2x^ sinS kP 
x 8(A- |fe + |)6(A- |fe_|) 



Q<,< 



A/s 



P 



w + ft, k + P cos 9 k p, k - —T7- sin 9 k p 



ipf I — oj, k — Pcosdkp, 7r + 6k 



P 
2Kl 



sin 9 



HP 



where f 1 = f°° ^ f 2 

Jk.uj J— 00 27r JO 

a factor of y/Kp into the definition of the fields, as is 
usual in the RG. The last two factors are to ensure that 
the actual momenta in the fermion lines are consistently 
within the cutoff. 



du r 27T d]h r' 

27T J -A 2-7T " 



rA dk 



We have absorbed 



where we have used the step functional* to set the cutoff 
to A/s. In this expression, we should now redefine {k 1 = 
sk, lu' = soj) to shift the cutoff back to its full value, and 
then use the definition of the field rescaling once again, 
to findi£ 



J 



s<,< 



= s 



-2 , 



I k' ,u/ 



■ sfi, k' + sP cos 9k' p, 9k 



Psin 9 k ' p ' 
2K F , 



\p[ i-u)\k' - sP cos 9 k - p. 



7T + 9 k ' + 



Psin 9 k ' p 
2K F 



This is the tree-level rescaling of the composite oper- 
ator; at one loop, the composite operator acquires an 
anomalous dimension, whose effects may be captured 
with another renormalization parameter Zq (we lump 
the factor of s~ 2 with this as well). Thus, we have the 
overall rescaling: 



Q<,< 



■A/2 



(22) 



where the ' denotes the fact that we have related fields 
in two different theories, the original and the one after 
the RG step. 

Following the same steps used in deriving the C-S 
equation previously, and combining the scaling of the 
descendant operator (Eq. I22[) with the inhomogeneous 
term (Eq. [2TJ) contributed by the pieces of the compos- 
ite operator lying above the cutoff, we have the following 
equation obeyed by the flow of the two-point function 
(OO): 



(Csl,pCfi,p)s,A 
1 

1 - 



27TWp 



(23) 

ZoZl(O s Q !S pO s n :S p)s s ,A(24) 



Once again, we may use the fact that the left hand 
side and thus the right hand side of this expression are 
independent of s to perform the same manipulations as 



before, and arrive at the expression 



1 



27n>F 
d_ 
ds 



2 + 2 7e > 



n,p/s,A 



O. 



o spCsn.sP 



(25) 



where we have defined 2 



We leave it as an exercise to the reader to show that 
the final term on the right hand side is equivalent to 
the action of the operator + P-§p on the two point 
correlation function; any discrepancies between the two 
vanish because of the symmetry of the Fermi surface. 
Making this substitution, we find 



n ii +p jp + ^T 9 +2+2 ^ 
1 



2ttvf 



P S,A 



(26) 



As a final step, we note that once again we need 
to remove a trivial delta function in going between 
the expectation value and the correlation function, 

(On,pOn.p) s ,A = To (to>Q) s n,n s P,-p- Thus ; we finally 
arrive at the inhomogeneous Callan-Symanzik equation 
for the two-point Cooper pair correlator, 



J 



rg- 21 (S!, 4 ; 9 ,A) = - 



27TWp 



(27) 
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IV. SOLUTION OF THE COOPER PAIR 
CALLAN-SYMANZIK EQUATION 

The next step is to solve the Callan-Symanzik equa- 
tions: given a bare coupling constant go, we wish to de- 
termine the long-wavelength, low-frequency behavior of 



,(0,2) 



Since the dimensionless combination of the fre- 



quency, momentum, and cutoff that enter the correla- 
tion functions must be of the form y/fl 2 + P 2 /A = fi/A, 
(where a phenomcnological velocity - that depends on the 
Fermi-liquid parameters such as m* - has been set equal 



to unityJi) Using this, we argue that + Pjjp- 



— A^- when acting on the correlation functions; this gives 
the slightly more tractable equation 
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r£' 2) (fi; 5 ,A) 



iTTVp 



(28) 



We may now solve this equation by the method of 
characteristics^, and find 



rg' 2) (fi;<?o,A) 



with x—g(x) 
ax 



-P(g(x)) and g(x = 0) = g 



1 



2-KVp J 1 x 



* dx 



^ e -2A x ^7o(9(y)) 



(29) 



Before we can complete our solution, we need to com- 
pute (3 and jq. We sketch the calculation in the Ap- 
pendix, and just cite the results here: 

P(g) = ~ag 2 
lo(g) = -ag (30) 

with a > 0. 

With these in hand, we hnd that 



We see that for go strictly zero, the expression diverges 
logarithmically as fl — > 0, reflecting the singularity in the 
zero-frequency, zero-momentum pairing response of the 
free Fermi gas at T = 0. However, for any finite 170, we 
find that the response is nondivergent: r^' 2 ' (Q; g, A) ~ 
l/2irvFgo as — > 0. Observe that, nicely enough, this 
answer itself diverges as g — > 0. 



go 



1 




2-nvpg 1 

where the flow of the coupling is given by 

H \ go 

ago log f 




(32) 



which clearly reflects the fact that g is marginal: as we 
take fl — > 0, g vanishes logarithmically. Note that we 
must always have Q < A, so that the logarithm in the 
denominator is positive and does not lead to any singu- 
larity as we take 17 — > 0. 

Since the fact that 9 f§J IS marginal, in the limit of 

interest, the second term in Eq. dominates, and we 
have 



rg' 2) (tt;s,A) 



1 


«l 


'ft 




2-Kvpgl 




)- 


\ a 








2ttv f 1 + ago 




.0 

' A 



go 



(33) 



V. ANOTHER ROUTE 

We have outlined our derivation above at some length 
for we were interested in a particular method of getting 
the answer, in which we follow the irrelevant coupling 
all the way to zero while continuing to renormalize. For 
the generic marginal coupling, this is the easiest way to 
go and even for fermions there are other problems, e.g. 
involving gauge fields, where we expect this technique 
will be the way to go. 

However, there is another route to our answer — as 
readers may guess by looking at it. In this approach we 
renormalize until we get to an exactly solvable problem 
and then we appeal to the exact results. In our prob- 
lem, the action with the purely BCS channel interaction 
corresponds to the reduced BCS Hamiltonian, which has 
infinite range interactions and is thereforw exactly solv- 
able by a saddle point computation in the infinite vol- 
ume limit. The same method shows that the RPA result 
for the superconducting susceptibility for this problem 
is exact. This has precisely the form (|3"3")l with A now 
being the scale at which we switch to the exact solution. 
Adding in the additional operator renormalizations gath- 
ered en route will change the answer but not the finiteness 
of the result or its behavior as go — s- 0. Indeed, in this 
approach it is also straightforward to explicitly include 
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the Landau couplings as the resulting Hamiltonian is still 
exactly solvable^. 

The general procedure we have described in this section 
is also what is used in the implementations of the the RG 
for interacting fermions known as the Functional Renor- 
malization Group (FRG), see for example Appendix B 
in Ref. [Tol The difference is that as in such work gener- 
ally relevant flows with multiple coupling constants get 
stopped at some scale the resulting problem is not typi- 
cally exactly solvable in a controlled sense. However, that 
has to do with the ends to which the FRG is put — the 
idea is the same. 



VI. DISCUSSION 



Our result (|33l) verifies the claim with which we 
began — namely that the superconducting susceptibility 
of the Fermi liquid is finite due to the intercession of 
the marginally irrelevant BCS coupling that is present 
even for repulsive interactions. This exact compensa- 
tion of the leading singularity by the irrelevant flow is 
somewhat surprising from the RG perspective, certainly 
if you compare with the results on ferromagnets in four 
dimensions that we reviewed in the introduction. Pos- 
sibly multi-band systems will lead to richer possibilities 
but that is a subject for future work. It is also of interest 
to extend the RG approach taken here to the derivation 
of the effects of the marginal flow on the electron Green's 
function. As this requires application of the RG at two 
loops, this will require going beyond the straightforward 
momentum shell method used in this paper. 
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Appendix A: Calculation of the Beta Function and 
Anomalous Dimension 



While we expect many of our readers to be familiar 
with the computation of the RG functions, for the sake of 
completeness, and because the anomalous dimension cal- 
culation may be unfamiliar to some, we briefly review the 



procedure here. Following Shankar—, we compute correc- 
tions to the action from integrating out modes between 
A — dA to A via the cumulant expansion. Essentially, 
we treat the coupling as a perturbation SS, and use the 
result 



a SS\ 



5 [<5S) + I((A-S) 2 )-((5S)) 2 



(Al) 



where all averages are over the modes being integrated 
out. In each term, we are to select the number of 'ex- 





FIG. 2: One-loop diagrams for the beta function and anoma- 
lous dimension of tptp 



ternal' legs, which will belong to the modes below the 
shell of integration, and this determines which term in 
the action will be corrected by the term obtained by in- 
tegrating out the remaining fields. In order to determine 
the anomalous dimension of the Cooper pair operator, 
the easiest method is to add a source term of the form 
Joty^ (and its complex conjugate) to the action and de- 
termine how it gets renormalized at one loop; a moment's 
thought will suffice to realize that this is equivalent to de- 
termining 7e>. 

The one-loop diagrams contributing to the /^-function 
and the anomalous dimension of O are shown in Fig. [5J 

2 dA 



Evaluating the first diagram, we find that 8g = ag 



'hi 

lis 



= —ag 



since 



or in other words that P(g) 

ds = — dA/ A . From the second diagram, we find that 
^j^- = ag^j-; combining this with the tree level rescaling 
we find that Jo — > s l ^ a9 Jo^ which gives us the result (cf. 
Eq. 122]) that 2 + 2^o = 2 — 2ag, or in other words, that 
70 = — ag. In these expressions, a is a positive constant, 
whose value is unimportant; the significant point is that 
it is the same constant in both (3 and 70. 
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